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Abstract 

Control and estimator gains are computed for linear-quadratic-Gaussian (LQG) op- 
timal control of the axial vibrations of a thermoelastic rod. The computations are based 
on a modal approximation of the partial differential equations representing the rod, and 
convergence of the approximations to the control and estimator gains is the main issue. 


1 Introduction 

The axial vibrations of a uniform rod are represented by a one-dimensional wave equation 
with constant coefficients, and thermoelastic damping in the rod is represented by a one- 
dimensional heat equation coupled to the wave equation. The solutions to the wave and 
heat equations are, respectively, the axial displacement and temperature fields in the rod. 
[See 1,2]. 

The length of the rod in this paper is normalized to 1. For active control, a single force 
is distributed parallel to the rod, uniformly over the portion sq < s < s\ of the rod. The 
equations of motion of the plant are then 


pwtt = (A + 2p)w S9 — a(3A + 2 p)6 9 + bu + br/i , £>0, 0<s<l, 


pc6t — k0 SB — 0qq:( 3A -(- 2/i)u>* 5 , t > 0, 0 < s < 1, 


where 


( 1, s 0 < s < s u 
1 0, otherwise. 
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(1.1) 

( 1 . 2 ) 

(1.3) 
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In these equations, w(t) = w(t,s) is the axial displacement, 0(t) = 9(t,s) is the temperature 
distribution and u(t) is the control force. We assume that the actuator force has the form 
u + ? 7 i where u(<) is the known control function and 771 is zero-mean Guassian white noise 
with intensity q\. The constants />, a, A,/i,# 05 c and k are physical constants with values to 
be given later. 

We assume that we have a sensor that measures the displacement at the left end of the 
rod segment over which the actuator force is distributed. This measurement is then 

y(t) = w(*,s 0 ) + % (1.4) 

where t)q is zero-mean Guassian white noise with intensity 1. 

In this paper, we use the boundary conditions 

w(t,0) = = 0, , _v 

0,(<,O) = 0,(M) = O, { ’ 

which mean that the rod is clamped and insulated at both ends. Because of the insu- 
lated, or Neumann, bondary conditions for the heat equation, zero is an eigenvalue of the 
open-loop system and the corresponding eigenvector represents a constant temperature dis- 
tribution. The span of this eigenvector is an uncontrollable and unobservable subspace, 
and the orthogonal complement of this subspace contains only states for which the average 
temperature along the rod is zero. Whatever the initial conditions and the control function, 
the average temperature in the rod therefore is a constant function of time. We will denote 
this average temperature by 8. 

2 The Control Problem 

We set 0 = 0 — 0 and define the state vector 

x(t) = (w(t),w t (t),6(t)). (2.1) 

We take the state space to be the Hilbert space 

E = 0, 1) x I 2 ( 0, 1) x I 2 ( 0, 1) (2.2) 

where Hq(0, 1) is the first-order Sobolev space containing functions that vanish at both ends 
of the interval. The system in (1.1) - (1.4) then has the form 

x = Ax + Bu + Btji, (23) 

y = Cx + Tfc, (2.4) 

where A generates a strongly continuous semigroup of contraction operators on E (see [3]). 
We note that |x(f)|!; is the sum of twice the mechanical energy in the rod and the integral 
over the rod of 0 2 at time t. 

The LQG optimal control problem in this paper is to find u to minimize 

J = lim S{L f ((Qx(t),x(t))^ + u 2 (t)]dt) (2.5) 

f— >oo t J0 
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where 


Qx = (w,u;t,0). ( 2 -6) 

The operator Q is chosen so that mechanical energy is penalized but temperature variations 
from 0 are not penalized. As usual, this problem separates into a deterministic linear- 
quadratic regulator problem and a state estimation, or filtering, problem. Each of these 
problems has a unique solution because the open-loop system can be shown to be uniformly 
exponentially stable. 

The optimal control law has the form 


u(t) = - Fx(t ) 

(2.7) 

where the state estimate x satisfies 


x = Ax + Bu + F(y — Cx). 

(2.8) 

The gain operators F and F are given by 


F = B*P 

(2.9) 

and 

F = PC* 

(2.10) 


where P and P are the unique nonnegative self-adjoint bounded linear operators on E 
satisfying the Riccati operator equations [4,5] 


A*P + PA — PBB*P + Q = 0 

(2.11) 

and 


AP + PA* - PC* CP + qiBB* = 0. 

(2.12) 

The gain operators can be represented in terms of elements of E\ i.e., 


FX = ((/l,/2)/3),z)£, (/l»/2,/3) . 

(2.13) 

F = (fi,hJ 3 )cE. 

(2.14) 


3 Approximation 

We approximate the infinite dimensional system in (2.3) and (2.4) with a sequence of finite 
dimensional control systems of the form 

in ~ A n X n "h (^*^) 

y n = C„x n + 7/o. (3-2) 

In [10], we compaxed two Galerkin approximations: a finite element scheme in which linear 
splines were the basis vectors and a modal scheme in which the open-loop eigenvectors of 
the distributed system were the basis vectors. The modal scheme gave faster convergence 
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for approximations to control gains like those in (2.9) and (2.13). We use the modal scheme 
here. 

It is easy to see that the cigenspaces of the open-loop thermoelastic rod with the bound- 
ary conditions in (1.5) are three-dimensional subspaces each spanned by a two-dimensional 
eigenspace of the undamped wave equation and a one-dimensional eigenspace of the heat 
equation. The eigenvectors are sine waves for the wave equation and cosine waves for the 
heat equation. Since the modal approximation amounts to projection onto a sequence of 
complete orthogonal subspaces, it is easy to show that the approximations to the open- 
loop semigroup and adjoint semigroup converge strongly, as commonly needed in numerical 
solution of infinite dimensional Riccati equations for distributed systems [4-7). 

For each n, we approximate the solutions to the infinite dimensional Riccati equations 
(2.11) and (2.12) by solving a pair of finite dimensional Riccati equations involving A n ,B n 
and C n . With the solutions to these matrix Riccati equations, we approximate the control 
and estimator gains as in [5]. In particular, we compute approximations to the functional 
gains (/i,/ 2 ,/ 3 ) and (/ 1 ,/ 2 ,/j) in (2.13) and (2.14). 

4 Numerical Results 

We chose the constants in (1.1) and (1.2) for an aluminum rod of length 100m (see [8, 9]). 
We normalized the length to 1, so that the constants take the numerical values 

p = 9.82 x 10~ 2 A = 2.064 x 10’ 1 p = 1.11 x 10 _1 
a = 1.29 X 10“ 3 4 c = 5.40 x 10" 1 
k = 7.02 x 10 -7 6 0 = 68. 

The actuator force is spread uniformly between so — .385 and s\ = .486, and the intensity 
of rji is qi = 1. 

Figures 1-3 show the approximations to the control functional gains (/i,/ 2 ,/ 3 ) for a 
range of approximation orders n, and Figures 4-6 show the approximations to the esti- 
mator functional gains (The number of modal subspaces used is n.) Because 

thermoelastic damping is very light in aluminum, as in most metals, many modes of the 
rod must be controlled actively. Therefore, many modal subspaces must be used in the 
approximations before the functional gains converge. 
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Figure 4(b); fl, the estimator gain, n « 18 - 33 
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